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Abstract 

The  first  Born  approximations  to  solutions  of  a  time- 
domain  integral  equation  were  used  to  obtain  the  backscat- 
tered  dilatation  wave  response  from  spherical  and  cylindri¬ 
cal  inclusions  of  arbitrary  homogeneous  anisotropic  elastic 
material  embedded  within  a  homogeneous  isotropic  host.  For 
large  ka,  where  the  validity  of  the  first  Born  approximation 
is  questioned,  the  time  waveform  responses  from  cylinders 
and  spheres  are  markedly  different?  for  ka  — >  0,  they  have 
identical  time  form  with  amplitudes  dependent  upon  the 
volume  of  the  scatterer.  Excellent  agreement  with  an  exper¬ 
imental  result  for  scattering  from  a  cylindrical  void  was 
obtained  for  a  value  of  ka  ■  0.32. 

A  "transparency  condition"  was  obtained,  allowing  that 
for  certain  combinations  of  both  density  and  stiffness  of 
the  scatterer  and  host,  the  scatterer  appears  transparent  to 
the  incoming  wave  in  the  first  Born  approximation. 

These  results  are  of  practical  significance 
destructive  inspection  of  fiber-reinforced  composite  materi¬ 
als,  with  elastic  waves  of  long  wavelength,  for  determining 
the  presence  of  porosity  remaining  in  the  composite  after 
manufacturing. 


vi 


Introduction 


I  . 


Background 

Fiber-reinforced  composite  materials  are  being  intro¬ 
duced  into  primary  structures  of  modern  jet  fighter  aircraft 
(e.g.  vertical  and  horizontal  stabilizers  on  the  F-15  and 
F-16  fighters).  These  materials  are  made  up  of  many  layers; 
each  layer  consists  of  many  rows  of  parallel  fibers  bound 
together  by  a  matrix  material  of  different  composition  than 
the  fibers.  Failure  often  occurs  in  these  composites  via 
delamination  which  initiates  at  gas  bubbles  (porosity) 
remaining  in  the  matrix  after  the  manuf acturing  process. 
For  this  and  other  reasons,  non-destructive  inspection  tech¬ 
niques  capable  of  determining  the  presence  of  the  porosity 
are  being  sought.  Current  efforts  include  inspection  with 
elastic  waves.  Items  of  interest  include  comparison  of  the 
response  from  the  fibers  to  that  from  porosity  and  the 
degree  to  which  the  fibers  may  appear  transparent  to  the 
wave  s . 

For  this  study,  the  fibers  are  modeled  as  cylinders  and 
the  porosity  as  isolated  spheres*  Many  studies  have  been 
done  on  the  scattering  from  cylinders  and  isolated  spheres 
[1-3],  linear  arrays  of  parallel  cylinders  [4],  and  two- 
dimensional  arrays  of  parallel  cylinders  [5-6] .  These  stu¬ 
dies,  however,  all  provide  frequency-domain  solutions. 
Current  inspection  techniques  utilize  broadband  time-domain 
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pulse  techniques;  thus,  a  time-domain  solution  for  the 
scattering  might  be  more  readily  applied  to  the  inspection 
problem . 

Problem  and  Scope 

The  problem  investigated  in  this  study  is  the  solution 
for  the  dilatation  wave  backscatter  from  cylindrical  and 
spherical  inclusions  resulting  from  incident  dilatation 
pulses  which  are  representative  of  the  pulses  produced  by 
piezoelectric  transducers.  Obtaining  and  exploring  the 
differences  between  backscatter  from  spheres  and  cylinders 
and  comparison  with  an  experimental  result  are  the  major 
objectives  of  this  thesis.  A  secondary  objective  is  a  study 
of  the  transparency  of  common  composite-reinforcing  fibers 
to  elastic  waves  in  order  to  assess  the  feasability  of  non¬ 
destructive  inspection  with  backscatter ed  dilatation  waves. 

The  analysis  is  limited  to  backscatter  from  single 
cylindrical  or  spherical  inclusions  of  arbitrary  homogeneous 
anisotropic  elastic  composition  embedded  in  a  homogeneous 
isotropic  host.  The  solution  for  the  cylinder  is  obtained 
for  normal  incidence  of  the  dilatation  wave  with  respect  to 
the  cylinder's  axis  of  symmetry.  All  media  are  considered 
to  obey  the  laws  of  linear  elasticity  and  are  assumed  to  be 
non-attenuative. 
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Approach  and  Presentation 

A  review  of  Lee's  development  [7]  of  a  time-domain 
integral  equation  for  the  scattered  field  is  given  in 
Chapter  II,  followed  by  a  statement  of  the  time-domain  first 
Born  approximation  and  the  form  of  the  solution  for  the 
scattered  pulse  resulting  from  this  approximation.  A  "tran¬ 
sparency  condition"  is  then  obtained  for  a  general  anisotro¬ 
pic  homogeneous  inclusion  in  a  homogeneous  isotropic  host. 
In  Chapter  III,  the  solution  for  backscatter  from  cylinders 
and  spheres  is  worked  out  in  detail.  Power  series,  Laplace 
transform  methods,  and  recursion  relations  are  used  to 
obtain  the  solutions.  Some  specific  results  in  the  form  of 
graphs  of  back scattered  responses  are  presented  and  dis¬ 
cussed  in  Chapter  IV,  where  comparison  to  an  experimental 
result  is  also  presented.  The  transparency  of  typical 
fibers  is  also  addressed  in  Chapter  IV.  Conclusions  and 
recommendations  are  presented  in  Chapter  V. 
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I I .  P 1 ane  Pulse  Scattering 

The  first  section  is  intended  to  be  a  brief  expose  of 
Lee's  [7]  theoretical  analysis  of  dilatation  wave  back- 
scatter  from  an  object  insonified  by  a  plane  dilatation 
pulse.  The  time-domain  first  Born  approximation  is 
presented  therein  to  linearize  the  resulting  expression  for 
the  scattered  displacement  field.  A  product  of  this  thesis 
is  presented  in  the  following  section  where  a  result 
obtained  in  (7]  for  scattering  from  a  void  is  generalized 
for  any  homogeneous  anisotropic  scatterer.  It  will  be  shown 
that,  in  the  first  Born  approximation,  certain  combinations 
of  density'  and  stiffness  of  the  host  and  scattering  materi¬ 
als  render  the  scatterer  "transparent"  to  the  incoming 
waves . 

The  first  Born  approximation 

Following  the  development  of  Lee  [7] ,  consider  a  piece 
R  of  linear  elastic  material  bounded  by  the  surface  dR .  A 
transducer  located  on  a  subset  A  of  9R  launches  a  dilatation 
pulse  which  encounters  a  region  B  of  different  linear  elas¬ 
tic  material  located  somewhere  within  R,  as  shown  in  Figure 
t.  If  B  is  far  enough  away  from  A,  and  if  the  radii  of  cur¬ 
vature  of  A  are  large  enough  compared  to  the  dominant 
wavelengths  of  the  pulse,  then  both  the  wave  incident  upon  B 
and  the  scattered  wave  received  by  the  transducer  can  be 
approximated  by  plane  waves.  The  analysis  which  follows  is 
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not  initially  restricted  to  baclocatter,  i.e.  the  observa¬ 
tion  vector  x.  may  point  in  any  direction.  (This  would,  of 
course,  require  another  transducer  to  receive  the  scattered 
waves.)  The  condition  for  backscatter  is  imposed  later  in 
the  chapter. 

The  material  densities  and  stiffness  tensors  are  given 
by 


f  %°  >  i  t  fc-B 

Y  $°  *  )  *  €-  B. 


(  1) 


and 


Qjkl^ 


CjjW  £  X  Sij  +  +  Sit  Sjk)  ;  X  ^  R-B 

C°ijVd  AC,-jW  ) 


,  \0  0  • 
where  A  and  ^  are  the  Lame  constants  in  R-B  and 

Kronecker  delta  defined  by 


the 


s 


l) 


)  J 
,  1  - ) 


(3) 


The  incident  pulse  is  of  the  form 


g  «X 


(4) 
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where  a  is  the  dilatation  wave  speed  in  the  region  R-B  and 
0 

e  is  a  unit  vector  characterizing  the  direction  of  propaga¬ 
tion  of  the  pulse.  The  constant  p  allows  the  time  origin  to 
be  adjusted  so  that  t  *  0  corresponds  to  the  instant  that 
the  pulse  first  encounters  B.  Thus, 


=  ~  VV\\IA 


u-*o 


(5) 


In  (4),  the  scalar  potential  $(s)  satisfies 


4(s^ 


e  C 


(6) 


and 


<j>ls)  =  0  ,  s  (7) 

A  minor  modification  of  equation  (5.4.2)  of  reference 
(7J  defines  the  pulse  length  Tp  as  that  time  after  which  a 
specified  fraction  £  of  the  total  energy  in  the  pulse 
remains  according  to 
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The  displacement  field  then  satisfies 


u'*c  (,x,0  s  -®ct<0 
+  ust  (.*  >-0  >  t>0 


(9) 


An  application  of  Love's  integral  identity  [8]  gives 
the  scattered  field  as 


ururi*  lip.* 

B 


l  4-bJ  dv; 


(  10) 


where  b^(x,t)  is  a  body  force  field  which  is  dependent  upon 
the  incident  displacement  field  along  with  the  density  and 
stiffness  perturbations 


(in 


In  (10)  the  dots  indicate  differentiation  with  respect  to 
time,  and  the  functional  u  is  given  by 
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(12) 


where 


r  ~  \  t\  -  I  i  ~y  I  (13) 

and  bQ  is  the  shear  wave  speed  in  R-B.  Equation  (10)  shows 
that  the  scattered  field  arises  in  an  obviously  nonlinear 
manner.  It  may  be  possible  to  solve  for  usc(x,t)  by  an 
iterative  solution  of  (10),  however,  a  simpler,  linear  prob¬ 
lem  results  by  assuming  that  the  interaction  of  the  scat¬ 
tered  field  and  its  derivatives  with  the  scatterer  is  much 
smaller  than  that  of  the  incident  field.  This  requires  that 

Acijkl'  J1*C  )  *  and  £  8C(x,t)  be  small  in  some  sense 

so  that 


9 


jSuJr,t|K,.XArM^] 

B 

B 

Equation  (14)  is  a  statement  of  the  "time-domain  first 
Born  approximation"  for  elastodynamic  scattering  which  then 
gives  the  scattered  field  as 

U\  ”  jjj  I  (15) 

& 

Useful  results  may  be  obtained  by  considering  far-field 
dilatation  wave  backscatter  resulting  from  a  plane  c  lata- 
tion  pulse  input.  Generally,  both  dilatation  and  shear 
waves  will  be  scattered  due  to  mode  conversion  at  the 
scatterer's  boundaries;  however,  at  distances  which  ere 
large  enough,  time-gating  may  be  employed  to  observe  only 
the  response  due  to  the  faster-traveling  dilatation  waves. 

By  expanding  equation  (IS)  with  (12),  discarding  terms 
which  are  0(x"2)  and  terms  corresponding  to  shear  waves,  and 
applying  Gauss'  divergence  theorem  to  what  is  left,  the 
far-field  scattered  dilatation  field  u*C(j,t)  is  obtained  as 
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where 


X  -  |  X 


(  17) 


and 


a  X 

X  =.  — —  (18) 

"  1*1 

Taylor's  expansion  of  r  and  order  of  magnitude  arguments 
allow  further  simplification  of  (16).  The  resulting  expres- 
sion  for  the  scattered  dilatation  field  for  an  incident 
plane  dilatation  pulse  is 


.  &  jj  j  ACKu,nW«e,S+^A^ 
0  ‘ 


i(< 


1  Qo  A* 


(  19) 
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where 


Y(V)  ==  i  C  t) 

(20) 


is  the  amplitude  of  the  incident  displacement  field,  and 


(21) 


A  homogeneous  anisotropic  scatterer 

The  result  obtained  in  equation  (19)  is  a  statement  of 
the  first  Born  approximation  for  the  most  general  anisotro¬ 
pic  inhomogeneous  scatterer  embedded  within  an  isotropic 
homogeneous  host.  Reference  [7]  goes  on  to  consider  further 
details  for  the  special  case  of  scattering  from  a  void.  In 
this  thesis,  equation  (19)  is  evaluated  for  a  homogeneous 
anisotropic  inclusion.  The  solution  has  the  same  form  as 
that  for  the  void,  with  the  addition  of  a  multiplicative 
amplitude  factor  which  depends  upon  the  materials'  elastic 
constants . 

The  perturbations  attributable  to  the  scatterer  are 
taken  as 

Z  o 

M  =•  -  $  (22) 
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and 


=  ”  C‘iM 


(23) 


B  B  0 

where  ^  and  the  0°^^  are  constants  within  B,  and  c^^i  is 

as  defined  in  equation  (2).  Substituting  these  into  (19), 

evaluating  with  the  Kronecker  deltas,  and  setting  x  *  -e  for 

backscatter,  yields 


(24) 


Consider  the  product  eie je^« ic i 1  in  (24).  If  the 
cijkl  ar®  contracted  according  to  the  method  described  by 
Nye  [9*131-149],  the  product  can  be  written  out  in  its  full 
glory  (for  future  reference)  as 


(25) 


(26) 


(27) 


Equations  (25)  and  (28)  together  give  the  backscattered 


dilatation  wave  amplitude  from  any  homogeneous  anisotropic 
scatterer  in  the  first  Born  approximation.  This  result  has 
the  same  form  as  that  obtained  for  a  void  [7] ,  i.e. ,  the 
backscattered  dilatation  wave  amplitude  is  proportional  to 

E> 

It  is  shown  in  (7]  that  performing  the  integration  results 
i  n 

TO  =  ^r  1  If  ait)  it 

L  L  (30a) 


or 


(30b) 


for  both  "|f(t)£  C  *  2  *  (  -oo ,  oo )  and  A(t)  £  2  *  ( -oo ,  oo )  ,  where 

equation  (30a)  is  obtained  for  -  0  and  (30b)  for 

» 

A ( 0 )  «  0  .  To  obtain  these  results,  a  coordinate  system  has 

been  introduced  with  the  z  axis  parallel  to  e  and  the  time  T 
is  defined  by 

r  - 

a0  (3i) 
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Note  that  )  -  0  so  that 

mi  n 


£  =  -2-*^;, 


(32) 


The  function  A (D  is  a  function  which  describes  the 
scatterer's  cross-sectional  area  as  a  function  of  twice  the 
one-way  transit  time  of  the  plane  wave  passing  through  the 
scatterer.  The  dependence  of  A(T)  upon  twice  the  travel 
time  makes  sense  physically,  since  the  backscattered  portion 
of  a  pulse  which  travels  some  distance  into  the  scatterer 
must  also  traverse  the  same  distance  back  through  the 
scatterer . 

By  combining  equations  (28)  and  (29),  the  scattered 
displacement  field  can  be  expressed  as 


— ' .. 
HIT  a*  x 


M 


(33) 


where  the  material-dependent  amplitude  M  is  defined  as 


n 


i  eKe,ewenCKUs 

gO  r  c° 

3  '• 


(34) 


It  is  instructive  to  consider  a  homogeneous  isotropic 
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In  terms  of  fami- 


inclusion  B  in  order  to  obtain  M. 

isot  ropic 

liar  engineering  quantities.  It  is  convenient  to  define 


=  r  +  a  a 

0 

A*  =  /*  +  Ar- 

=  e°  f  a£ 

v*  =  +-  av 


<  35) 


(  36  ) 


(37) 


(38) 


where  E  and  V  are  Young's  modulus  and  Poisson's  ratio, 
respectively.  Substitution  in  (34)  of  an  equation  for 

g 

Cijkl'  *n*l*9ous  to  equation  (2)  for  isotropic  media,  leads 
to 


it 


f  ,  C* 

+  7* 


-1 


(  39) 


Note  that  this  expression  reduces  to  the  result  obtained  in 
(7]  for  scattering  from  a  void,  namely 


MvoU  -  -  2-  ,40’ 

Equation  (39)  is  alternately  expressed  in  terms  of  E  and  V 

as 
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M  =  £_  ,  )0-^° )  _  O 

1  '\so4ro^c  (|-D°  )(.'  *-VB)(|-lV*  ) 


(41  ) 


which  is  linear  in  Ae  and  A))  . 


Consider  the  ratio  in  (41) 


Substituting  for  VB  from  (38)  and  expanding  yields 


R.  = 


2Av> 

|-xv>* 


If  A V  satisfies 


M  «  \ 


ltV° 


-  -v* 

7.  v 


-1  <,»°  <~'A 

-  '/^  <V°  <  '4 


(42) 


(43) 


(44) 


then  Rv  can  be  approximated  by  keeping  the  first  two  terms 
in  the  geometric  series  for  the  last  two  factors  in  (43) 
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(45) 


1  + 


2.P*  {l-v°) 


Thus,  M  is  also  linear  in  Ay  for  small  . 

The  possibility  of  a  scatterer  appearing  "transparent" 
to  a  plane  dilatation  pulse  in  the  first  Born  approximation 
is  suggested  by  equations  (34)  and  (39).  By  setting  M  *  0, 
one  form  of  the  "transparency  condition"  for  an  isotropic 
scatterer  is  obtained  from  (41)  as 


2~ 


(46) 
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Ill . 


The  results  of  the  previous  chapter  for  a  scatterer  of 
arbitrary  shape  are  applied  to  spheres  and  cylinders. 
First,  the  cross-sectional  areas  of  the  scatterers  as  a 
function  of  time  are  determined.  A  model  for  the  amplitude 
of  an  incident  displacement  pulse  as  a  function  of  time  is 
then  introduced.  The  form  chosen  for  the  pulse  is  not  only 
mathematically  attractive,  it  also  closely  models  the  pulse 
which  is  produced  by  physical  transducers.  With  the  area  and 
pulse  functions  determined,  <^SC(t)  is  evaluated  according  to 
equation  (30a)  for  spherical  and  cylindrical  inclusions  of 
arbitrary  isotropic  linear  elastic  constants.  The  analytic 
character  of  these  solutions  is  investigated,  followed  by  a 
look  at  scattering  in  the  long  wavelength  limit. 


Cross-sectional  scattering  areas 

Consider  a  plane  dilatation  wave  pulse  traveling  in  the 
positive  *  direction  and  incident  upon  a  spherical  scatterer 
B,  of  radius  s,  as  shown  in  Figure  2a.  The  cross-sectional 

A 

area  A  (z)  "seen"  by  the  wave  is  the  area  of  the  family  of 

3 

circles  subtended  by  the  plane  on  the  sphere 


r ,  -snss 

|  o  ,  o^U«.»*vJise  (47) 


Recalling  equations  (5)  and  (31)  for  p  and  f,  the  area  Ag(^) 
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f 

Figure  2.  Plane  wave  interaction  with  a  sphere 
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is  obtained  as 


Ag(r)  is  parabolic  as  illustrated  in  Figure  2b. 

Now,  consider  a  dilatation  wave  pulse  incident  upon  a 
cylinder,  where  the  vector  e  is  perpendicular  to  the  axis  of 
symmetry  of  the  cylinder,  as  shown  in  Figure  3a.  If  the 
radius  of  the  cylinder  is  c,  and  the  length  of  the  cylinder 
over  which  the  wave  can  be  approximated  by  a  plane  wave  is 
L,  then  the  cross-sectional  area  is  the  family  of  rectangles 
described  by 

-C  i  3  i  C 

'  (49) 


oi’Ui*’  ut  i? 


(50) 


which  describes  half  of  an  ellipse,  as  illustrated  in  Figure 
3b. 

Since  <j>ac(t)  is  given  in  equation  (30a)  as  a  convolu¬ 
tion  of  a  pulse  amplitude  function  with  M<£>),  Figures  2b  and 
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3b  suggest  (by  coniidtrina  graphical  convolution)  that, 
perhaps,  spheres  and  cylinders  will  give  similar  scattered 
time  waveforms.  This  topic  is  considered  in  the  section  on 
long  wavelength  scattering.  However,  the  finite  slope  of 

the  sphere’s  area  function  at  r  ■  0  and  T  m  —  compared  to 

*0 

the  infinite  slope  of  the  cylinder's  area  function  at  the 
corresponding  points  will  be  shown  to  give  characteristi¬ 
cally  different  scattered  responses  when  the  long  wavelength 
limit  is  not  valid.  Note  that  it  is  precisely  the  infinite 
slope  A£(  0 )  which  requires  the  use  of  equation  (30a)  instead 
of  (30b)  to  obtain  $sc(t)  for  the  cylinder. 

It  will  prove  convenient,  in  what  follows,  to  introduce 
a  normalized  time 


t  ’ 


(5  1) 


where 


b 


,  sphere 
,  cylinder 


(52) 


This  allows  the  incident  pulse  to  be  characterized  in  terms 
of  the  width  of  the  scatterers'  cross-sectional  area  time 
functions.  Note  that  2b  is  the  time  it  takes  the  wave  to 
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completely  traverse  the  scatterer  in  the  incident  direction 
and  again  in  the  backscattered  direction. 


Using  (51)  and  (52)  in  (48)  and  (50),  the  normalized 
area  functions  are 

f(0Lo  ^  (t'-t'*)  n  Oi  t‘  i  I 

n  oVktruiJC 

W  (53) 


Astf)  = 


and 


AS?)  = 


ZWLa,-/  t  -T,u  ,  0  -i’i 


O  VWt'  u  i  s « 


(54) 


Incident  pulse  model 


A  mathematically  tractable  and  physically  representa¬ 
tive  model  of  the  incident  plane  dilatation  pulse  is  given 
as 


0  , 


(55) 


where  u)  is  the  center  frequency  of  the  pulse.  The  factor  of 
t  is  used  to  obtain  ^(0)  -  0  for  use  of  equation  (30a). 
Differentiating  twice  with  respect  to  time 
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r 


-2.^^  ***wt  t  co^cit'  ,  i:  >o 


YW* 


0 


>  t*0 


(56) 


Using  (51)  and  (52)  in  (55)  and  (56),  the  normalized 
incident  pulse  amplitude  functions  are 


YlfcV 


bt '  $fv»  O't'  t*  >  0 


(57) 


and 


2\>(^~u){  o‘i’  - e^coWf 

~2|{  ^  sUw‘4  •v-'los^2^  CGw't'  ^  t*>c 


Y«V 


,t$0 


(58) 
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where 


CO'  -  2 bcO 


(  59) 


is  the  normalized  radian  frequency. 


T 


In  terms  of  the  period 


T  ' 


(60) 


gives  the  period  of  the  pulse  relative  to  the  time  width  of 
the  scatterers '  area  functions. 

It  is  revealing  to  put  this  time  domain  pulse  normali¬ 
zation  in  the  perspective  of  the  frequency  domain  wavenumber 
normalization  ka  employed  throughout  the  literature,  where  a 
is  the  characteristic  dimension  of  the  scatterer,  and  k  « 

ao 

is  the  wavenumber.  If  a  corresponds  to  either  radius,  s  or 
c,  and  recalling  (52)  and  (60)  for  b  and  T*, 


ka 


(61) 


Solution  for  the  sphere 


with  all  the  pieces  now  available,  the  normalized  solu* 

O' 

tion  for  the  backscattered  amplitude  function  <p  ® c  ( t  *  )  for 
scattering  from  a  sphere  may  be  obtained.  Either  form  of 

equation  (30)  may  be  used  to  obtain  $*c(t’)  since  ( 0 )  *  0. 
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For  the  form  given  by  (30b),  the  required  second  derivative 
of  the  area  function  is 


TTaH1  [Str'}  -  1  -  ,  o<t'£| 


The  delta  functions  greatly  simplify  the  convolution  with 
the  incident  field  in  (30b)  and  only  the  convolution  with 
the  constant  (-2)  must  be  worked  out  in  detail.  The  situa¬ 
tion  is  not  so  pleasant  for  the  cylinder  and  equation  (30a) 

is  required  due  to  the  singular  behaviour  of  A  (t*)  at  T1  »  0 

c 

and  at  f  «  1.  In  order  to  provide  a  unified  treatment  of 
scattering  from  both  objects,  the  form  of  (30a)  is  chosen, 
obtaining  the  solution  in  terms  of  the  area  function  con¬ 
volved  with  second  derivative  of  the  incident  pulse. 


Substituting  equations  (53)  and  (58)  into  (30a)  pro¬ 


duces 


xiz'-r'jjr  ,t'r  0 


where  dr  “  2bdC  is  employed,  and 


oC 


> 


O  <t'<[ 
t'  >  \ 


1 


1 


(  64) 


which  limits  the  integration  to  that  portion  of  the  pulse 
which  is  encountered  by  the  wave  at  time  t'.  Introducing 
the  complex  exponent 


=  tJoq.  t  i  W* 


(65) 


and  using  Euler’s  formulae  for  the  sine  and  cosine  func¬ 
tions,  equation  (62)  is  expanded  and  rearranged  to  obtain 


'  fcv  I  - 
re  6t 


+i 

-t  UlO 


,»  *V 


-')  A  10)  r  e.  dr' 
u 


i  to 


(66) 
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where 


AlO* 


and 


+  ifcV 


-<+t' 


(67) 


BIO* 


i  i 


e 


(68) 


The  asterisk  denotes  complex  conjugation  and  i  =  \/-T .  The 

integrals  in  (66)  are  easily  handled  by  defining 


CjM  s  *)  Z'J  ^ 1  dr’ 


(69) 


Integration  of  (69)  by  parts  leads  to  the  recursion  relation 


CjW  = 


j  K.V 
•<  e. 


_  (\  (*\ 

\c+  J‘l 


(  70) 


which  is  particularly  useful,  since  C^fat)  is  readily 
evaluated.  Note  also  that 


C  i  ^  t  j  i 

C\W  =  )  t  e  dr 


(71) 


so  that  (66)  can  be  rewritten  as 
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=  <  [/ICOC,W-A‘U,Krw]  t 

-I  [gU'KM-BVK’w] 

+1  [BU'KiW-fcVK’W] 

+  <  [Atf)CjW-^(^C‘W]  (  i‘>o 

(  72) 

If  Im( • )  denotes  the  imaginary  part  of  a  quantity,  the  back- 
scattered  wave  amplitude  from  a  sphere  is  finally  obtained 
as 


2TT  a,4  \ y1' 


-2  |l,[Aft')c.w]+’ 

4>  CiC«i)j 


(73) 
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where,  recall,  «,  equals  t'  or  1  according  to  equation  (64). 

While  equation  (73)  is  not  particularly  revealing,  it 
is  a  beautifully  simple  result  that  is,  together  with  (69) 
and  (70),  quite  amenable  to  numerical  evaluation.  It  is 
easy  enough  to  observe  from  (67)  through  (70)  and  (73)  that 
<j>®c(t')  has  the  same  center  frequency  as  2f(t')  and  falls  off 
as  t'e  for  large  t'. 


So lution  for  the  cylinder 


By  a  similar  development  as  that  for  the  sphere, 

<psc(t')  is  obtained  for  the  cylinder  by  substituting  equa- 
*  c 

tions  (54)  and  (58)  into  (30a)  producing,  upon  rearrangement 


-A^cosoi'*1  +  v>o)’f  y 

■  ^  VT'-rl  S"n  o'  r'dt’ 


go’£'  -  Z-k'^u)  OStJ'i 

*  (  e  t' VT^r71  cosu’r'  dr' 

Jo 

^  t'1  -1  ‘i-is  <^<*)  5  ia  u  t ' 


x(  sf  z''X,x  S'nut'  t’>0 

;o  ’ 


^d1 
>e  ® 


(74) 


where,  again  o<  ■  t '  or  1  according  to  equation  (64)  .  The 

integrals  in  (74)  lead  to  a  solution  for  ^®c(t’  )  which  is 

/•» 

not  as  straightforward  as  <j)®c(t')  obtained  for  the  sphere. 
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Consider  the  integral 


CoStoV  <Jx‘ 


(75) 


Expansion  of  the  radical  in  a  binomial  series  gives 

X,iO  =  \ 

o 

provided  0<t'<1  .  Now,  since  the  series  is  uniformly  con¬ 

vergent  within  the  interval,  the  summation  and  integration 
may  be  interchanged  to  obtain 


(77) 


where 


« 

Z  ®  CoscJ  r  <iT 


Simi larly , 


5 

)o 


(78) 


(79) 


3 4 


become  s 


where 


-it 


lT’ 


S'*  o  T  dr 


(80) 


(  81  ) 


The  functions  Ac(t'  )  and_/\.®(t‘)  lend  themselves  to  recursion 
3  J 

relations  which  make  them  computationally  attractive. 
Integrating  (78)  and  (81)  by  parts  yields 


^  ''A.+) 


1*0* 


+  u  L 


f^cosw 


r 


+•  w  $i*i 


c St' 


_  1 _ as  (.,•) 


(82a) 


and 


1'^  r  T 

-  I  2.W*  S\*i  tO't  -  cj'coSwH' 

(aJb^V+-u'1  L  b  J 

A5  (*’)  _  -iii _ a*  «• ) 


(82b) 


The  appearance  of  the  recursion  index  j  in  the  numerators 
of  the  forward  recursion  relations  (82)  is  not  desirable 


i 


since  truncation  errors  introduced  in  the  numerical  evalua¬ 
tion  of  the  functions  will  be  amplified  as  j  gets  large.  A 
numerically  stable  evaluation  scheme  results  from  manipula¬ 
tion  of  (82)  which  leads  to  the  reverse  recursion  relations 


AV)  =  -j-L—  [  i'  (.O' 


(83a) 


-,  +  1  L 


(83b) 


Note  that 


limA^(t*)  -  limA®(t’)  *  0  (84) 

j-»oo  3  j-»oo  3 

which  is  easily  seen  by  considering  (78)  and  (81)  remember¬ 
ing  that  0 <t  '  <  1 . 

Because  of  (84),  the  reverse  recursion  relations  are 
also  convenient  in  that  an  exact  expression  is  not  needed 
for  one  of  the  recurring  functions  in  order  to  start  the 
recursion.  For  large  j. 


A3_,  >  j-*oo 


(85) 


I 


36 


i 


37 


which  can  be  written  as 


2 

Now,  let  v 


ku 


This  gives 


£ 

<ik 


Furthermore,  the  substitution  v  *  ifj,  where 
renders 


which  is  recognized  to  contain  a  form  [10:297]  of 
plex  error  function 


X« 


£ 


-i  J&erf(iSkT) 


Performing  the  indicated  differentiation  [10:298] 


Back  to  eva  luating  ( t '  )  and/lgt*’),  Euler’s 


(89) 


(90) 


i  -  Z=T, 


(9  1) 


the  com- 


(92) 


(93) 


f ormu lae 
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for  the  sine  and  cosine  functions  in  (78)  and  (81),  and  an 


application  of  (93)  to  both  yield 


(94b) 

where,  as  before  K—  -  2bg  +  iu)’.  Equations  (74),  (78),  and 
(81)  can  now  be  combined  to  give  the  solution  valid  for 
0<t*<1,  i.e.  when  the  front  edge  of  the  pulse  is  within  the 
scatterer.  The  result  is 
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(95) 


where 


A, a')  = 

-  U  J Swut*  Coso’V 

(96a) 

AtuO  = 

£  ^  ]  C.0  Su  ‘  V+jlk  ^  t  -  Uu  j  Si  v\  cj' t 

(96b) 

AjU'): 

Smg o'i’  -  co^o’i' 

(96c  ) 

Cesci’t'  -V-  2.^^  SMU't1 

(96d) 

The  f  unct  tone  /\*  ( t  •  )  end  J\  *(t'i 


ere  sterted  by  equetion 
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(86),  evaluated  by  equation  (83),  and  compared  to  (94)  to 
insure  a  desired  degree  of  accuracy. 

The  solution  for  t'^1  must  still  be  obtained.  Equation 

(74)  for  4sc(t')  is  equivalent  to  a  form  more  indicative  of 
1  c 

the  convolution  operation  (indicated  by  an  asterisk),  and 

thus,  amenable  to  Laplace  transform  methods.  The  desired 

form  of  £  sc ( t '  )  is 
T  c 


■fLCfc 


(97) 


where 


-  (t'e^Slv.u't'^Xr  9(t‘)  <98a) 

(98b) 

0Ct’)  (98c) 


^Ct')  =  COS  Jt')  * ©W ) 


(98d) 


I 
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and  where 


r 


ecn  - 


Ji'-i"  ,  Olk'i.  I 

O  \  oV^rutie. 


(99) 


If  M*1  denotes  the  one-sided  Laplace  transform  of  a  quan¬ 
tity,  the  required  transforms  are  evaluated  in  terms  of  the 

transform  variable  s  as 


L[te^COwt 


L[ 


]  = 


Co  Suit  j  — 


and  f  rom  (11:138) 


to 


[^'wl 


U) 


V  w 


J±i_ 


\eVv  I,W 


kZ 


( 100a) 


(  100b) 


( 100c) 


(  lOOd) 


(101) 


where  I^(.)  is  the  modified  Bessel  function  of  the  first 
kind  of  order  V  •  Denoting  L  [  f  ( t  *  )  ]  by  F^s),  for  i  « 
1,2,3  ,4,  the  F  ^  (  s  )  are  obtained  by  the  multiplication  of  the 
appropriate  transforms  in  (100)  by  the  transform  in  (101). 
The  f  (t*)  are  conveniently  obtained  by  summing  the  residues 

st  • 

at  the  poles  of  e  F^(s),  which  poles  are,  in  all  cases. 


(102) 


Note  that  - - -  is  an  entire  function  of  s  and,  therefore, 

s  *  0  is  not  a  pole  of  the  F^(s).  Saving  the  details,  the 
result,  valid  for  t*>1,  is  obtained  as: 


Inn  = 

4LO> 


in[m 


+  [Ett1)] 

-  t*  I„  [x>U')] 


(103) 


where 
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(104) 


and 


+u'-i\  X.(£ 


(105) 

and  where  Re (  «  )  denotes  the  real  part  of  the  quantity. 

Equations  (95)  and  (103)  together  give  the  solution  for 
scattering  from  a  cylinder  for  all  t',  including  the  criti¬ 
cal  point  where  the  solutions  are  pieced  together  at  t*  •  1. 
Again,  the  form  of  the  solutions  is  not  particularly  reveal¬ 
ing,  however,  the  same  observations  as  those  made  for  the 
scattering  from  the  sphere  can  be  made  with  respect  to  the 
scattered  center  frequency  and  the  form  of  the  solution  for 
large  t ' . 

Given  the  simplicity  of  the  solution  for  t’>1  compared 
to  the  solution  for  (Kt*<1,  one  is  tempted  to  look  for  a 
Laplace  transform  solution  valid  for  all  t'.  However,  the 
above  inversions  are  valid  only  for  t*>1.  Consider,  for 
example,  e8t  F3(s)  in  the  form 

U’p.  utr  I.  (VO 

e  '5  Z  M06) 

The  asymptotic  expansion  of  Iv(z>  {10t377]  for  large  z  and 
|arg(z)|  along  with  the  fact  that  I^z)  is  an  odd  function 
of  z  yields 
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-*/». 

—  e. 


M— *oo  4- Tr  s 


(107) 


Thus , 


,4-'  con-  e. 

est  p-.U'i  »  — — —  - - - 

3  is.i-*oo  i_}-irs  s  (.s-sp4Us-s>f'> 


)  2-  <cu^U}<3—  (108) 


S  t  * 

Since  the  poles  of  e  P^(s)  are  in  the  left  half  of  the 

g  t  * 

complex  plane,  e  F3 ( s )  is  required  to  vanish  as  is|->oo 

in  order  to  sum  residues  to  obtain  f  (t*).  However,  for  s  in 

the  left  half  plane  (with  a  negative  real  part), 
s  t  * 

lim  e  F  (s)  *  oo,  unless  t*>1,  which  therefore  limits 
I  s  I  -  >  oo 

the  Laplace  transform  solution  to  the  region  t*>1. 


Smoothness  of  the  solutions 

The  interaction  of  the  incident  pulse  with  scatterers 
having  abrupt  onset  and  terminating  boundaries  raises  ques¬ 
tions  about  the  smoothness  of  the  scattered  waveforms. 
Recall,  the  scattered  field's  second  time  derivative  must  be 
sufficiently  integrable  so  that  equation  (14)  is  satisfied 
for  the  first  Born  approximation  to  be  valid.  The  solution 
for  the  cylinder  is  particularly  suspect  due  to  the  infinite 
slope  of  the  area  function  Ac (r' )  presented  to  the  wave  at 

C'  ■  0  and  t*  *  1*  Thus,  4>8C(t*)  is  treated  here  first,  in 

c 

the  form  presented  by  equations  (97)  and  (98). 

By  writing  out  the  convolutions  as  integrals,  keeping 
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* 


the  first  terms  in  the  expansions  for  all  factors  in  the 
integrands,  and  then  performing  the  integrations,  the  result 
obtained  is 


Thus , 


■  Vi 


1  «•  i 


(109) 


H  U,H 


t  cct 


(110) 


%  t‘"‘A 

2_ 


(111) 


Note  that  <t)sc(t,),  5sc(f  )  ,  and  (j>8c(t’)  are  all  zero  if 
T  c  T  C  T  c 

•« 

t  ’  <  0 .  Thus,  only  (j5®c(t’)  is  discontinuous  at  the  onset  of 

-1/2 

scattering,  with  an  infinite  discontinuity  from  at 
s ingu la  r ity • 


Recall,  in  the  convolutions  of  (97),  the  upper  limit  of 
integration  is  t’  for  0<t * < 1  .  Differentiating  with  respect 


i 
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to  t*  according  to  Leibnitz'  rule  [12]  yields 


+  ^iijt'-e3-  ,o£V$i 

(112) 

where  the  f^(t')  obviously  give  back  combinations  of  the 
fi(t')  and  the  last  term  results  from  differentiating  f^(t') 
ala  Leibnitz.  For  t'j>1» 

-^±1  =-zio^\(-^-Uw']>UV 

+•  J\>  4-  logoi'  ^  U') 

(113) 

where  the  g^Jt’)  are  the  same  convolution  integrals  as  the 
f^(t’),  except  that  the  upper  limit  of  integration  is  1, 
instead  of  t' .  Since  gt<1)  -  f^l),  for  i  ■  1,2, 3, 4,  com- 

parison  of  (112)  and  <113)  implies  chat  |'c<t’l  is 


♦a“tf  \  _ 
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continuous  at  t*  «  1.  Moreover,  since  the 

f^(t' ) £  C*  (-00,00) ,  the  important,  physically  expected, 

result  that,  <4sc(t')£  C *  1  * ( -o© , o© )  is  obtained.  Now  dif- 
1  c 

ferentiating  <^*c(t')  introduces  a  factor  of  ( t 1 -t  *  2  )  ~  1//2  for 

•• 

0<t'<1  which  obviously  means  that  CD  3  c  ( t  *  )  is  not  continuous 

1  c 

at  t '  *  1 .  In  fact , 


1  im 
t '  1 


4bLa2 

o 


-oo 


(114) 


Similarly,  for  t">1 


lim 


t  ’  -*  1 


+ 


4bL*o 


<  oo 


(115) 


Thus,  ^®c(t' )  has  an  infinite  discontinuity  at  t*  »  1  just 
as  it  does  at  t'  -  0. 

It  is  reasonable  to  question  the  validity  of  the  first 

Born  approximation  in  the  light  of  these  results  for 
•• 

^^C(t'),  since  this  relates  directly  back,  through  equation 

(33),  to  equation  (14).  The  approximation  is  not 

•  1/2 

threatened,  however,  since  the  t'  singularity  is  suffi¬ 

ciently  integrable  in  (12)  and  volume-integrable  to  insure 
that  ( 14 )  holds . 
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The  same  analysis  applied  to  the  sphere  yields  that 


3T  Wu)t'  ,  t'<-Ct 


Again,  the  s 
of  scatteri 
tinuity  b<*i  i 


econd  derivative  is  discontinuous  at 
ng  although  finite  for  the  sphere, 
s  alternately  expressed  as 


the  onset 
The  discon- 


t> 


cO 


2LWcl 


(117) 


where  ka  is  as  mentioned  in  the  section  on  the  incident 
pulse  characteristics.  The  second  derivative  is  found  to  be 
discontinuous  at  t’  *  1  by  -bW,  and  more  importantly, 
^®c(t'  )  £  C( n) ( -00,00) . 

Long  wavelength  1 lmlt 

The  time-domain  first  Born  approximation  can  already  be 
thought  of  as  a  long  wavelength  approximation  in  its  own 
right.  This  seems  reasonable  since  it’s  an  approximation 
based  upon  weak  interaction  of  the  incident  field  with  the 
scatterer;  weak  interaction  is  intuitively  appealing  for 
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wavelengths  long  compared  to  the  size  of  the  scatterer. 
However,  it  appears  that  the  approximation  as  stated  in  (14) 


could  be  satisfied 


by  sufficiently  small  perturbations 


and  Ac,  ,  even  in  shorter  (or  at  least  resonant)  wavelength 
i]kl 

regimes.  This  certainly  is  an  item  begging  for  further 


attention . 


The  long  wavelength  limit  is  stated  according  to  T’>>1, 
or  equivalently,  ka<<^\  In  this  case,  equation  (30a)  is  par¬ 
ticularly  revealing  since  both  A^Ct)  and  A c(t)  now  sample 

the  second  derivative  of  the  incident  pulse  over  suffi¬ 
ciently  small  intervals  so  as  to  produce  results  similar  to 
convolutions  with  a  delta  function.  Thus,  in  the  long 
wavelength  regime,  spheres  and  cylinders  produce  essentially 
identical  backscattered  fields  which  have  pulse  amplitudes 
given,  very  nearly,  by  the  second  derivative  of  the  input 
pulse  amplitude.  This  result  is  obtained  independent  of  the 
shape  of  any  small  scattering  object.  The  most  significant 
difference  between  the  backscattered  fields,  in  the  long 
wavelength  limit,  will  be  the  amplitudes  of  the  scattered 
fields  due  to  the  "strength"  of  the  delta  function.  This 
strength  depends  upon  the  size  of  the  object  and  is  mani¬ 
fested,  in  the  case  of  cylinders  and  spheres,  in  the  volumes 
2  3 

L/Tc  and  4f(s  /3 .  Of  course,  the  material  dependent  ampli¬ 
tude  factor  M  introduced  earlier  plays  a  role  in  the  scat¬ 
tered  amplitude. 

As  an  example  of  how  the  delta  function  strength  and  M 
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affect  a  backscattered  response,  consider  the  problem  of 
distinguishing  the  response  from  a  lonq  cylindrical  inclu¬ 
sion  and  a  spherical  void  in  the  lonq  wavelenqth  limit.  The 
scattered  time  functions  will  have  indistinguishably  similar 
form,  to  within  a  minus  sign  attributable  to  M,  and  only  the 
relative  amplitudes  may  be  strikingly  different.  The 
cylindrical  inclusion  may  very  well  have  an  M  which  is  quite 
small  compared  to  that  for  the  spherical  void.  However,  if 
the  length  L  over  which  the  cylinder  is  illuminated  by  the 
approximately  planar  wave  is  large  enough,  the  response  from 
the  weakly  scattering  inclusion  can  dwarf  the  response  from 

the  relatively  strongly  scattering  void,  due  to  the  strength 
2 

LTfc  which  is  unbounded  in  L.  This  problem  is  of  consider¬ 
able  interest  to  the  Air  Force  for  non-destructive  inspec¬ 
tion  of  fiber-reinforced  composite  materials,  where  detec¬ 
tion  of  gas  bubbles  (porosity)  in  the  epoxy  matrix  (host)  is 
desirable.  It  would  be  nice  if  the  cylindrical  fibers  were 
relatively  transparent  to  the  incident  dilatation  pulse,  as 
compared  to  a  spherical  void,  allowing  the  response  from  the 
porosity  to  stand  out  against  a  "background"  of  cylinders' 
response.  Further  attention  is  given  this  item  in  the  fol¬ 
lowing  chapter. 
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IV.  Numerical  Results 


Some  of  the  numerical  techniques  used  for  computer 
evaluation  of  the  scattered  responses  from  cylindrical  and 
spherical  inclusions  are  presented.  Responses  obtained  for 
incident  pulses  having  various  normalized  center  frequencies 
ka  are  then  presented  to  illustrate  the  differences  and 
similarities  between  scattering  from  cylinders  and  spheres. 
An  experimental  result  for  scattering  from  a  cylinder  is 
presented  and  a  comparison  made  to  the  response  predicted  by 
the  theory.  Finally,  the  transparency  of  the  fibers  in  a 
fiber-reinforced  composite  material  is  addressed. 

Numerical  techniques 

Recall,  the  incident  pulse  amplitude  is  modeled  as  an 
exponentially  damped  sinusoid  according  to  equation  (55). 
For  the  following  results,  the  length  of  the  pulse 
defined  by  equation  (8)  was  chosen  to  be  the  time  after 
which  ten  percent  of  the  energy  in  the  pulse  remains.  With 
a  desired  pulse  length  and  center  frequency  specified,  sub¬ 
stitution  of  (55)  into  (8)  leads  to  a  transcendental  equa¬ 
tion  for  the  damping  q  which  was  solved,  most  con¬ 
veniently,  by  the  simple  method  of  bisection,  as  described 
in  [13:65]  . 

The  error  function  and  modified  Bessel  functions  of 
complex  arguments  were  evaluated  by  routines  based  upon  the 
forms  of  the  functions  given  in  [10] .  Power  series  and 
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continued  fractions  were  utilized  to  provide  at  least  eight 
significant  figures.  Asymptotic  expansions  were  employed  in 
the  appropriate  regions.  The  reverse  recursion  relations 
(83)  were  considered  to  provide  accurate  results  if  they 
matched  the  closed  form  expressions  (94)  to  seven  figures. 
Finally,  251  points  were  calculated  for  each  theoretical 
plot  displayed  from  here  on. 

Time-domain  backscatter  waveforms 

A  typical  incident  pulse  and  its  second  derivative  are 
shown  in  Figure  4,  illustrating  the  particular  case  when  the 
pulse  length  is  equal  to  the  period.  This  may,  of  course, 
be  specified  for  a  pulse  of  any  period.  The  specific  case  of 
scattering  when  T*  -  0.5  (ka  «  if)  is  shown  in  Figure  5.  The 
most  interesting  difference  between  the  responses  from  the 
cylinder  and  sphere  occur  around  the  point  t '  ■  1 .  It  is  at 
this  time  that  the  portion  of  the  wave  which  makes  the  round 
trip  through  the  entire  scatterer  is  felt  at  the  observation 
point  where  the  wave  first  enters  the  scatterer.  The  more 
interestingly  "structured"  response  from  the  cylinder  is  due 
to  the  more  abrupt  nature  of  the  change  in  density  and 
stiffness  which  the  wave  sees  as  a  result  of  the  infinite 
slopes  of  the  cylinder's  area  function.  While  the  first 
Born  approximation  is  not  guaranteed  (or  expected)  to  be 
valid  when  ka  is  this  large,  these  results  are  interesting 
in  that  they  predict  marked  differences  in  the  backscatter 
from  cylinders  and  spheres  at  shorter  wavelengths.  Figure  6 
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b. 

Figure  4.  Incident  pulse  waveforms 
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illustrates  the  scattering  when  ka  *  ft,  for  the  same  pulse 

shape  (  T  m  t*  )  as  that  in  Figure  4.  These  results  illus- 
P 

trate  less  structure  around  t'  -  1  than  for  ka  ■  TT.  This  is 

to  be  expected  since  increasing  T  and  T*  have  the  effect  of 

P 

narrowing  the  bandwidth  of  the  frequency  spectrum  of  the 
incident  pulse  and  centering  it  around  a  lower  center  fre¬ 
quency.  This  effect  then  carries  through  to  the  spectrum  of 
the  scattered  pulse  and  is  evidenced  in  a  time-domain  pulse 
of  less  structure. 

Comparison  to  experimental  result 

While  the  first  Born  approximation  may  be  suspect  for 
values  of  ka  as  large  as  in  the  examples  above,  excellent 
agreement  between  an  experimental  result  and  theoretical 
predictions  for  smaller  ka  will  be  presented. 

The  experiment  conducted  at  the  Air  Force  Materials 
Laboratory  is  illustrated  in  Figure  7.  A  voltage  waveform 
f(t),  input  to  the  transducer,  produces  a  plane  dilatation 
pulse  which  interacts  with  a  cylindrical  void  (radius  50.8 
t*m  )  located  inside  a  piece  of  linear,  homogeneous,  and  iso¬ 
tropic  elastic  material.  The  scattered  wave  gives  rise  to 

an  output  voltage  g  (t).  The  incident  wave  also  reflects 

sc 

from  the  planar  "back  wall"  producing  the  response  g.  (t). 

DW 

This  response  can  be  considered  to  be  unperturbed  by  the 
incident  and  reflected  waves  traveling  through  the 
scatterer,  since  the  magnitude  of  the  scattered  field  is 
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Figure  7.  Scattering  Experiment 


much  smaller  than  the  incident  and  reflected  fields*  It  is 
shown  in  [14]  that  treating 

jscW  “  W  <118) 

is  completely  equivalent  to  treating  the  convolution  (30a) 
for  the  scattered  field. 

Figure  8  shows  agreement  between  theory  and  experiment 
for  scattering  from  the  cylindrical  void  when  ka  *  0.32 
(T*  »  4.85).  The  solid  curve  in  Figure  8a  is  the  back  wall 
reflection,  modeled  numerically  by  the  broken  curve  given  by 

7(t)  ~  (119) 

The  scattered  response  and  the  numerical  prediction  are 
similarly  illustrated  in  Figure  8b.  Such  agreement  between 
theory  and  experiment  lends  credence  to  this  approach  to 
elastic  wave  scattering.  The  lower  center  frequency  of  the 
reflected  pulse  can  be  attributed  to  high  frequency  attenua¬ 
tion  present  in  the  host  material  and  not  accounted  for  in 
the  model.  The  scattered  pulse  is  not  affected  as  much, 
since  it  travels  through  only  about  one-fourth  the  distance 
that  the  reflected  pulse  travels  through.  Figure  8c  shows 
the  second  time  derivative  of  equation  (119).  Figures  b  and 
c  illustrate  the  result  predicted  earlier  that  the  scattered 
field,  in  the  long  wavelength  regime,  is  very  nearly  given 
by  the  second  time  derivative  of  the  incident  field's  pulse 
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Figure  8.  Comparison  of  an  experimental  result  (solid  line)  with  the 
theoretical  prediction  (broken  line)  for  backscatter  from  a  cylinder 
for  ka  ■  0.32.  The  length  of  the  abscissae  is  1  psec  for  all  graphs. 


amplitude . 


Trans parency  considerations 


The  transparency  of  a  single  carbon  or  graphite  fiber 
in  an  epoxy  matrix  is  considered  based  upon  the  fiber  elas¬ 
tic  constants  reported  by  Smith  [15] .  Adopting  the  notation 
in  [15] ,  the  coordinate  axes  are  defined  with  the  3  axis 
parallel  to  the  fiber.  For  a  plane  dilatation  pulse  nor¬ 
mally  incident  upon  the  fiber,  e3  «  0.  Characteristics  of 
the  fiber  (modeled  as  having  hexagonal  crystal  symmetry) 
imply  that  c  «  c2g  -  0  ,  cn  -  c22  ,  and 

c„,  *  (c«.  -  c  „  _  )  .  Using  these  to  evaluate  (25)  for  the 
66  11  12 

product  ^e^e^c^ ,  the  amplitude  Mfiber  is  obtained 

according  to  (34)  as 


M 


(120) 


Note  that  this  is  the  same  result  obtained  for  isotropic 
media  in  (39).  This  is  a  consequence  of  the  hexagonal  model 
of  the  fiber  and  the  normal  incidence  of  the  dilatation 
pulse.  It  is  well  known  [16:116]  that  a  hexagonal  structure 
will  appear  the  same  to  a  normally  incident  wave,  indepen¬ 
dent  of  where  _e_  lies  in  the  1-2  plane,  i.e.,  the  fiber  is 
"transversely  isotropic." 

Table  1  illustrates  the  scattering  amplitudes  for  the 
fibers  considered  in  [15]  . 
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Table  I 


Fiber  Constants  and  Scattering  Amplitudes 


Fiber 

f  3 
(g/cm  ) 

cn 

(  1010pa) 

Mf iber 

H 

WYB 

1.32 

3.13 

5.10 

6.50 

T-  2  5 

1.38 

1 .94 

2.89 

2.09 

T-  40 

1  .57 

1  .42 

2.06 

1.06 

T-50 

1.67 

1.25 

1.82 

0.83 

T-50S 

1  .69 

1  .22 

1.78 

0  .79 

T-75S 

1  .88 

0 .97 

1.47 

0 .54 

VYB 

1.53 

4.92 

8  .69 

18.88 

PAN 

1  .72 

2.31 

3 .88 

3.76 

HTS 

1  .69 

1.60 

2  .5  1 

1 .58 

T-400 

1  .77 

- -  . 

2.3  3 

3  .97 

3.94 

The  epoxy  density  was  taken  as  -  1.16  g/cm3  and  c® 1 
was  calculated  from  E°  *  3.27x109  Pa  and  V°  “  0.35  according 
to 

_  E  (j-y) 
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The  last  column  in  the  table  compares  the  backscattered 
energy  from  the  fibrous  material  to  that  from  a  void  of  the 
same  shape.  For  only  three  materials  is  the  backscattered 
energy  less  than  for  a  void. 

Practical  consideration  of  detecting  porosity  in  an 

array  of  fibers  must  also  include  the  geometry  of  the 

scatterer,  i.e.,  the  "strength"  of  the  scatterer  mentioned 

in  the  section  on  long  wavelength  scattering.  The  factors 
3  2 

4ft s  /3  and  flLc  are  what  is  to  be  considered,  not  the  fac- 
2  3  2 

tors  4bLaQ  and  2iraQb  appearing  in  equations  (73),  (95),  and 

(103),  which  are  due  to  the  time  normalization  employed. 
Consider,  for  example,  fibers  8  microns  in  diameter, 
illuminated  by  a  normally  incident  wave  which  is  approxi¬ 
mately  planar  over  a  length  of  1/4  inch  (the  size  of  a  typi¬ 
cal  transducer  face).  Without  considering  the  material 
parameters,  the  amplitude  response  from  a  single  sphere  will 
exceed  that  of  a  single  fiber  for  porosity  diameter  greater 
than  about  3.3  mils  (1  mil  »  25.4  microns  ).  By  including 
the  material-dependent  amplitudes  from  Table  1  and  requiring 

that  (4*Vs3/3)M  be  greater  than  Lc2M. . .  ,  one  finds 

void  tlDer 

that  the  lower  limit  of  porosity  size  that  produces  a  larger 
backscattered  response  than  an  8  micron  diameter  fiber 
varies  little  compared  to  the  range  of  energy  ratios  in  the 
table.  For  energy  ratios  varying  from  0.54  to  18.88,  the 
lower  limit  of  "detectable"  porosity  varies  from  3  to  5.5 
mil  in  diameter. 
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V .  Conclusions  and  Recommendations 

Conclusions 

The  first  Born  approximations  to  solutions  of  a  time- 
domain  integral  equation  were  used  to  obtain  the  backscat- 
tered  dilatation  wave  response  from  spherical  and  cylindri¬ 
cal  inclusions  of  arbitrary  homogeneous  anisotropic  elastic 
media  embedded  within  a  homogeneous  isotropic  host.  For 
relatively  large  values  of  normalized  center  frequency  ka, 
the  time  waveform  responses  from  cylinders  and  spheres  are 
markedly  different  due  to  the  difference  in  slopes  of  the 
scatterers '  cross-sectional  areas  which  the  wave  sees  as  it 
first  meets  the  scatterer  and  then  as  it  exits  the 
scatterer.  For  ka  — »  0  the  responses  from  the  cylinder  and 
sphere  have  identical  time  form  and  are  given  by  the  second 
derivative  of  the  input  pulse  time  profile;  they  differ  in 
amplitude  and  the  difference  depends  upon  the  volume  of  the 
scatterer.  Excellent  agreement  with  an  experimental  result 
for  scattering  from  a  cylindrical  void  was  obtained  at  a 
value  of  ka  »  0.32. 

A  "transparency  condition"  was  obtained,  which  states 
that  for  certain  combinations  of  both  density  and  stiffness 
of  the  scatterer  and  host,  the  scatterer  appears  transparent 
to  the  incoming  wave  in  the  first  Born  approximation.  For 
inclusion-host  combinations  which  do  not  satisfy  the  condi¬ 
tion  exactly,  a  useful  quantity  is  the  energy  scattered  from 
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an  inclusion  normalized  to  that  scattered  from  a  void.  This 
is  of  practical  significance  for  non-destructive  inspection 
of  fiber-reinforced  composite  materials  with  long 
wavelengths.  A  limitation  exists  on  the  size  of  spherical 
voids  which  give  a  larger  amplitude  response  than  a  cylindr¬ 
ical  inclusion  of  a  given  size  based  upon  both  the  relative 
transparency  of  the  inclusion  and  the  volume  of  both 
scatterers  illuminated  by  the  wave.  For  an  Bp  diameter 
fiber  illuminated  over  1/4  inch,  this  lower  limit  on  the 
spherical  void  size  ranges  from  3  to  5.5  mil  in  diameter  for 
a  wide  variety  of  fiber  compositions.  Since  these  values 
are  within  the  range  (1-20  mil)  of  practical  interest,  it 
appears  from  this  simple  analysis,  that  an  inspection  tech¬ 
nique  based  solely  upon  ba ck scattered  ^wave s  should  not  be 
dismissed  without  further  analysis. 

Recommendations 

Extensions  of  this  study,  immediately  applicable  to  the 
non-destructive  inspection  problem,  should  include  con¬ 
siderations  of  backscatter  from  spheres  in  the  proximity  of 
an  infinite  linear  array  of  parallel  cylinders.  All  angles 
of  incidence  with  respect  to  the  plane  normal  and  the  direc¬ 
tion  parallel  to  the  fibers  should  be  considered,  since  nor¬ 
mal  incidence  alone  is  not  usually  used  in  current  inspec¬ 
tion  techniques.  Investigation  of  the  effects  of  periodic, 
almost  periodic,  and  random  spacing  of  the  array  could  prove 
useful.  Since  the  number  of  cylinders  needed  to  simulate 
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the  problem  numerically  would  be  very  large  in  the  long 
wavelength  regime  of  interest,  the  computer  program 
developed  for  use  in  this  thesis  would  not  be  useful,  pro¬ 
viding  as  it  does,  complete  solutions  for  the  scattering 
from  each  fiber  in  the  array.  Further  analytical  work,  sum¬ 
ming  the  responses  from  the  cylinders,  should  be  carried  out 
first,  and  the  scattering  from  spheres  obtained  against  a 
"background"  response  from  the  cylinders.  The  extension  to 
many  randomly  distributed  spheres  near  the  array  would  then 
provide  a  model  which  closely  approximates  the  physical  sys¬ 
tem  in  the  non-destructive  inspection  problem. 

Basic  questions  regarding  the  validity  of  the  time- 
domain  first  Born  approximation  need  to  be  answered.  The 
statement  of  the  approximation,  equation  (14),  suggests  that 
perhaps,  larger  ka  problems  might  be  solved  by  this  tech¬ 
nique  provided  the  perturbations  in  density  and  stiffness 
presented  to  the  wave  are  sufficiently  small.  Certainly,  as 
Aj>  — *  0  and  — >  0,  the  backscatter  goes  to  zero,  and  a 
regime  of  first  Born  validity  should  exist  for  small 
material  perturbations.  Another  nagging  point  is  the  artif¬ 
icial  unboundedness  of  the  material-dependent  amplitude  fac¬ 
tor  M  in  the  ratios  of  density  and  stiffness.  An  upper 
bound  needs  to  be  obtained  for  the  magnitude  of  M  beyond 
which  the  first  Born  approximation  ceases  to  be  usable. 
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